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Abstract. We consider computational and implementation issues for the comple- 
tion of monomial sets to involution using different involutive divisions. Every of 
these divisions produces its own completion procedure. For the polynomial case 
it yields an involutive basis which is a special form of a Grobner basis, generally 
redundant. We also compare our Mathematica implementation of Janet division to 
an implementation in C. 



1 Introduction and Basic Definitions 

In our previous paper [|J we described our first results on implement- 
ing in Mathematica 3.0 || different involutive divisions introduced 
in HJi|,||]; the completion of monomial sets to involution for those di- 
visions and application to constructing Hilbert functions and Hilbert 
polynomials for monomial ideals. 

In the present paper we pay more attention to efficient computation 
and propose some algorithmic improvements. Besides, we shortly de- 
scribe an implementation of Janet division in C and compare the run- 
ning times for both implementations. Though in this paper we consider 
involutivity of monomial ideals, all the underlying operations with in- 
volutive divisions and monomials enter in more general completion 
procedures for polynomial HQ and differential systems ||. 

Let N be a set of non-negative integers, and M = {x^ 1 • ■ ■ xf^ \ di £ 
N} be a set of monomials in the polynomial ring K[x\, . . . ,x n ] over 
a field K of characteristic zero. By deg(u) and degi(u) we denote the 
total degree of u E M and the degree of variable X{ in u, respectively. 



For the least common multiple of two monomials u, v G M we shall use 
the conventional notation lcm(u,v). An admissible monomial ordering 
is denoted by y, and throughout this paper we shall assume that it is 
compatible with 

x\ y x 2 y ■ ■ ■ y x n . (l) 

Definition 1. [|| An involutive division L on M is given, if for any 
finite monomial set U C M and for any u G U there is given a sub- 
monoid L(u, U) of M satisfying the conditions: 

(a) If w G L(u, U) and v\w, then v G L(u, U). 

(b) If u, v G U and uL(u, U) n u L(« , U) + 0, 
then « G vL(v, U) or w G uL(u, U). 

(c) If v G f7 and u G {/), then f7) C L(n, C7). 

(d) If y C [/, then L(u, U) C L(u, V) for all mGY. 

Elements of L(n, C/) are called multiplicative for u. If w € uL(u, U) we 
shall write and call u an (L-)involutive divisor of tu. In such an 
event the monomial v = w/u is multiplicative for n and the equality 
w = uv will be written as io = n x v. If n is a conventional divisor of 
io but not an involutive one we shall write, as usual, w = u ■ v. Then 
v is said to be nonmultiplicative for u. 

For every monomial u G U, Definition [l] provides the separation 

{ Xl ,...,x n } = M L (u, U) U NM L (u, U), (2) 

Ml(u, U) n NMl(u, U) = 0, of the set of variables into two subsets: 
multiplicative Ml(u, U) C L(u, U) and nonmultiplicative NMl(u, U)n 
L(u, U) = 0. Conversely, if for any finite set U C M and any u G U the 
separation (||) is given such that the corresponding submonoid L(u, U) 
of monomials in variables in Ml(u, U) satisfies the conditions (b)-(d), 
then the partition generates an involutive division. 

Definition 2. || Given an involutive division L, a monomial set U is 
involutive with respect to L or L— involutive if 

(Vn G U) (Vw G M) (3v G U) [ uw G 17) ] . 

In this paper as well as in |H] we shall consider the following eight 
different involutive divisions studied in ||[||5|] : 



Example 3. Thomas division ||. Given a finite set U C M, the variable 
Xi is considered as multiplicative for u EU \idegi{u) = max{degi(v) \ v G 
U}, and nonmultiplicative, otherwise. 

Example 4- Janet division Q]. Let the set U C M be finite. For each 
1 < i < n divide U into groups labeled by non-negative integers 
di,...,ck' 

[di, . . . , di] = { u eU | dj = degj(u), 1 < j <i }. 

A variable Xj is multiplicative for u 6 XJ if i = 1 and deg\{u) = 
max{degi(v) \ v £ U}, or if i > 1, u E [d% , . . . , c4_i] and degi{u) = 
max{degi(v) \ !)£ [d x , . . . , 

Example 5. Pommaret division ||] . For a monomial u = xf 1 ■ ■ ■ x^ fe 
with dk > the variables Xj,j > k are considered as multiplicative 
and the other variables as nonmultiplicative. For it = 1 all the variables 
are multiplicative. 

Example 6. Division I ||. Let U be a finite monomial set. The variable 
Xi is nonmultiplicative for u 6 XJ if there is t> G C/ such that 

x il ■ ■ ■ x< iZ u = lcm ( u i v )i 1 <m< [n/2], dj > (1 < j < m) , 
and Xi G . . . ,x im }. 

Example 7. Division II Q. For monomial u = x^ 1 • • • xf n the variable 
Xi is multiplicative if di = d max (u) where d max {u) = max{d\, . . . , d n }. 



Example 8. Induced division |Q. Given an admissible monomial order- 
ing >- a variable Xi is nonmultiplicative for u G U if there is v G £7 
such that v -< u and degi{u) < degi(v). 

To distinguish these divisions we use the abbreviations T, J, P, I,II,D 
In the implementation described below, three orderings are used to 
induce division in Example ||: lexicographical, degree-lexicographical 
and degree-reverse-lexicographical. For these three induced divisions 
we shall use the subscripts L, DL, DRL, respectively. 



Every of the above divisions generates its own procedure for com- 
pletion of a monomial set to involution by means of its enlargement 
with involutively irreducible nonmultiplicative prolongations. Given a 
monomial basis and an involutive division, the following algorithm 
MinimallnvolutiveMonomialBasis || produces the uniquely de- 
fined minimal involutive basis of the ideal. 



Algorithm MinimallnvolutiveMonomialBasis: 

Input: U, a finite monomial set 
Output: U, the minimal involutive basis of Id(U) 

begin 1 

U := Autor educe(U) 2 

choose any admissible monomial ordering -< 3 

while exist u G U and x € NMl(u, U) s.t. 4 

u ■ x has no involutive divisors in U do 5 

choose such u, x with the lowest u ■ x w.r.t. -< 6 

U:=UU{u-x} 7 

end 8 

end 9 



Here Autor educe(U) stands for the conventional (non-involutive) au- 
toreduction. 

2 Implementation Issues 

In this section we will describe some observations that allow to speed 
up the steps of the algorithm MinimallnvolutiveMonomialBasis 
significantly. Some of them are applicable to different divisions, oth- 
ers are concerned with the completion procedure in general. The basic 
operations on monomial sets are the same for the computation of in- 
volutive bases of polynomial HQ and differential systems ||, so the 
improvements described here are relevant for these computations. 

Our package provides a framework for studying the effect of us- 
ing different divisions and optimizations. It is implemented using a 
"generic programming" approach which allows to start with a straight 
forward implementation of the algorithm and introduce more efficient 
procedures for special situations later. 

The following statement returns the minimal involutive basis of a 
monomial set U with respect to Janet division and with lexicographic 
selection ordering: 



minimallnvolutiveMonomialBasis [Janet] [U, lexOrder] 



To extend the package for a new involutive division (called, say, 
newDivision), one would only have to write the specific version of the 
function separation which computes the multiplicative and nonmul- 
tiplicative variables of a monomial u £ U w.r.t. the set U: 

separation [newDivision] [u_, U_] : = . . . 

All the other steps in the algorithm would then be executed by func- 
tions that are generically defined for any involutive division. 

On the other hand, an optimized procedure for a specific situation 
can be introduced later to override the generic version. The pattern 
matching mechanism in Mathematica dispatches to the specific version 
wherever it is appropriate. 

Monomials are represented as multiindices, i.e. the monomial x % i ■ 
. . . ■ x l ™ is represented as the list of its exponents . . . ,i n }. Thus, 
the set U = {u%, . . . ,u m } can be considered as a m x n— matrix of 
integers. For every monomial u, we use two additional lists of length 
n: a list giving the separation of the variables for u, and a similar list 
containing notes about the prolongations that have already been done. 

We will now describe observations that can be used to make the 
basic operations of the algorithm MinimallnvolutiveMonomialBa- 
sis faster. Functions like 1cm will be applied also to multiindices, with 
the obvious meaning. The set notation is used for lists, assuming that 
the order of the elements is given somehow. U = {u\, . . . , u n } is a list 
of monomials, and u is always an element of U. 

The first step is to compute the separation for each of the input 
monomials. For globally defined divisions, this is done irrespective of 
the other monomials in U. For Janet division (Example we made 
use of the following remark: 

Remark 1. When the list U is sorted lexicographically in decreasing 
order, the groups [di, . . . , di\ mentioned in the definition are grouped 
together. These groups are sorted lexicographically with respect to 
their labels of any fixed length i. The sorted list starts with the group 
labeled [d\ ma x], d\ max = max degiU, the monomials in [d\ max ] have x\ 
as a multiplicative variable. We can split the list into groups given by 



labels of length 1 and proceed recursively within each of them, next 
considering degrees in the second variable X2, and so on. 

For a division Dy (Example ||) that is induced by some ordering 
>-, we can use an auxiliary list: 

Remark 2. Let the monomials be sorted in descending order: u\ y 
. . . y u n . We call the elements of the list cm(U) := {mi, . . . , m n \rrn = 
lcm(uj, . . . , u n ),i = n, . . . , 1} the cumulated multiples of U. By defi- 
nition, variable Xj is nonmultiplicative for U{ if and only if it has a 
higher degree in m^: degjUi < deg^rrii. Thus, all we have to do is com- 
pute the list cm(t7) of cumulated multiples and then compare each 
u £ U against its corresponding entry in cm(U). 

For Division I, we are not aware of any property that would allow 
us to accelerate the computation of separations in a manner similar to 
Janet or Induced divisions. 

The following observation can be used to speed up the process 
of finding a minimal nonmultiplicative prolongation (line 6 of the al- 
gorithm). Let us denote the minimal (w.r.t. the chosen ordering >-) 
nonmultiplicative prolongation by a given variable x with Py(x). 

Remark 3. Let U be sorted w.r.t. the completion ordering: u\ y . . . >- 
u n . Let Uiandx be fixed such that Uj • x is a minimal nonmulti- 
plicative prolongation w.r.t. >-. Then m ■ x is an element of the set 

{Py( Xl ),...,Py(x n )} . 

This follows directly from the minimality of Uj • x. Furthermore, 
Ui is the minimal monomial having X clS cl nonmultiplicative variable, 
because v ■ x y u ■ x implies v y u. 

The remark obviously extends to the more general situation of the algo- 
rithm, where some of the nonmultiplicative prolongations have already 
been considered. We keep a list P = {Py(xi), . . . , Py(x n )} of nonmul- 
tiplicative prolongations, one for each variable x\, . . . , x n , sorted by 
the completion ordering. Let v = Ui ■ Xj be the minimal prolongation. 
It is removed from P and checked for involutive divisors. If v is in- 
volutively reducible, we have to add another prolongation w.r.t. the 
same variable Xj to P. Otherwise, we add v to the monomial set and 
recompute the separations and P. 

The next step in the algorithm is to search for an involutive divisor 
w of a nonmultiplicative prolongation v = u-x. In the polynomial case, 
the efficiency of this search can be even more important, since we may 



want to involutively reduce every term of a prolonged polynomial. 
Recall that for an involutively reduced set U, there can be at most one 
such w. We present now some optimizations that apply to increasingly 
specialized situations. 

The following remark uses a special property of involutive divisions, 
taking into account that v is a nonmultiplicative prolongation of an 
element of U. 

Remark 4- Let U be an involutively autoreduced set of monomials and 
v = u-x a nonmultiplicative prolongation of some u G U. If a monomial 
w G U is an involutive divisor of v then deg x w = deg x v. 

Since u ■ x should be involutively reducible by w, we can write 
u ■ x = w x (u ■ x/w). If w = v = u ■ x, we are done. If w ^ u ■ x and 
w\u, then u = w x (u/w), which contradicts our assumption that U is 
involutively autoreduced. 

One can gain even more by considering particular divisions. Consider 
a Janet-autoreduced set U. Let us denote the longest common prefix 
of two monomials u, v by lcp(u, v), where lcp(u, v) := (ui, . . . , u^) with 
(ui, . . . , Ufc) = (vi,..., Vk), and k the maximal index for which and 
Vk coincide. If u\ ^ v%, we define lcp(u,v) := (). More generally, we 
use lcp(v, U) to denote the longest common prefix that v shares with 
some monomial from the set U. 

Remark 5. Assume that we search for a Janet - involutive divisor w of 
a monomial v. Then, w is in the class C defined by the label lcp(v, U). 
Let lcp(v, U) = (vi, . . . , Vk)- Every involutive divisor w = (w\, . . . , w n ) 
is also a conventional divisor, thus wi < Vi,i = l,...,k. We show 
by contradiction that Wi = Vi for i = 1, . . . , k. Let s be the smallest 
integer 1 < s < k such that w s < v s . Then, x s is nonmultiplicative 
for w because there exists a monomial in the class (v\, . . . , which 
has higher degree in x s , and w is not an involutive divisor of v. 

Note that this remark applies to arbitrary monomials v, not only those 
resulting from a nonmultiplicative prolongation. 

Consider a nonmultiplicative prolongation v = u-x. For Pommaret 
division, an involutive divisor w is reverse lexicographically greater 
than u. For a division that is induced by y, either u-x = woruyw 
holds. 

These properties together with Remark 3.12 in |jj suggest that one 
should keep the monomials sorted with respect to some suitable order, 
and use this order as completion order, too. 



Finally, when we find no involutive divisor, we have to add the 
prolongation to the set and adjust separations for all monomials ac- 
cordingly. Except for globally defined divisions, this step is potentially 
very time consuming. 

Remark 6. For all divisions discussed so far, the following holds for a 
monomial u G U: NM(u, U U {v}) = NM(u, U) U NM(it, {u, v}). 

A detailed discussion of this fact can be found in ||]. After adding a 
monomial v to U, this remark allows us to compute only the "pairwise" 
separations for every u G U. 

Specific divisions give rise to more improvements. 

Remark 7. Let v be a monomial, and assume that v has no involu- 
tive divisor in the Janet-autoreduced set U. Then, the separation may 
only change for monomials in the class lcp(v, U) = {v\, . . . , The 
separation of the variables x\,. . . ,Xk is left unchanged. Furthermore, 
the separation of the variables x±, . . . , x\. for the new monomial v can 
be copied from the separation of any of the monomials in the class 
lcp(v, U). 

Remark 8. Consider adding a nonmultiplicative prolongation v = u-Xj 
to an autoreduced set w.r.t. some induced division Dy- 

Only the variable Xj can change from multiplicative to nonmul- 
tiplicative, and it can do so only for monomials s >- v satisfying 
degjS = degjV — 1. 

Not all of the improvements mentioned here were actually implemented 
in the package. Our experience suggests that sometimes the practical 
performance in Mathematica differs from what one expects from look- 
ing at the algorithm. This is due to the interpreted nature of Mathe- 
matica and its flexible evaluation mechanism. Operations which are 
performed in the kernel are usually much faster than their equivalent 
expressed in a user defined function, and it was often a matter of trial 
and error to decide which variant of an operation one should use for a 
given division. 

In practice, the size of the resulting involutive basis is certainly the 
dominating factor for the overall running time of the algorithm. It was 
thus worthwhile to invest more programming work in improvements 
for those divisions which yield relatively small involutive bases (see 
below). 



The improvements for Janet division resulted the biggest gain in 
speed compared to the generic implementation. When the completion 
ordering is lexicographic, all optimizations described above are applied. 
For induced divisions Dy , we always use >- as completion ordering and 
Remark ^ to recompute the separations. Only for Division I, the time 
for changing the separations dominates the time for the other basic 
operations. Division I is also the only division for which the property 
mentioned in Remark |B] is used. The optimizations for finding an in- 
volutive divisor described above have a positive effect for all divisions. 

We have applied the package to examples taken from various sources. 
For each polynomial system, we computed the degree - reverse - lexico- 
graphical Grobner basis and took the resulting set of leading monomi- 
als as input to the algorithm MinimallnvolutiveMonomialBasis. 
As we described in Q the output can then be used to compute the 
Hilbert function, the Hilbert polynomial and the index of regularity of 
the corresponding polynomial ideal. 

Example 9. Consider anxn matrix A = (ctij) n n with unspecified 
entries. The condition A 2 = leads to a system of n 2 polynomial 
equations in the variables an, . . . , ai n , a2i, . . . , a nn . We treated the 
leading monomials of the degree reverse lexicographic Grobner basis, 
where the variables are ordered according to an >-...>- a\ n y a 2 \ >- 

• • • r Ot nn . 

Example 10. The system of "n -th cyclic roots" is a well known exam- 
ple. For n = 4, it is given by: 

xi + x 2 + x 3 + x 4 = 0, 

X\X 2 + X 2 X 3 + X 3 X 4 + X4X1 = 0, 

xix 2 x 3 + X2X3X4 + X3X4X1 + X4X1X2 = 0, 
X1X2X3X4 — 1 = 0. 



Example 11. The Reimer system in 5 variables: 

1 - 2x\ + 2x\ + 2x 2 3 + 2 *x 2 4 - 2x\ = °> 
1-2x1 + 2x 2 + 2xl + 2*x\ - 2x1 = 0j 
1 - 2x\ + 2x\ + 2x\ + 2 * x\ - 2x\ = 0, 



1-2x1 + 2x \ + 2x1 + 2* x\ - 2x\ = Oj 
l-2x\ + 2x\ + 2x\ + 2 * x\ - 2x% = 0. 



Example 12. The Katsura system in 7 variables: 

x\-x x + 2x\ + 2x1 + 2x \ + 2x l + 2x l + 2x 7, 

2x 2 xi + 2x 2 x 3 + 2x3X4 + 2x4X5 + 2x 5 x 6 + 2x 6 X7 - x 2 , 

2x3X1 + 2x2X4 + 2x3X5 + 2x4X6 + 2x5X7 + x\ — X3, 

2x4x1 + 2x2x5 + 2x3X6 + 2x4x7 + 2x2x3 — X4, 

2x 5 xi + 2x 2 x 6 + 2x3X7 + 2x 2 x 4 + x\ — x 5 , 

2x 6 xi + 2x 2 x 7 + 2x 2 x 5 + 2x3X4 - x 6 , 

xi + 2x 2 + 2x3 + 2x4 + 2x5 + 2x6 + 2x7 — !■ 



The following table shows the results of applying the algorithm Mini- 
mallnvolutiveMonomialBasis to our examples. In the first three 
columns, the size of the input is given where m is the number of 
monomials, n is the number of variables, and d is the maximum to- 
tal degree of the input monomials. The divisions are indicated by the 
abbreviations used above. For each division, we give the length of the 
minimal involutive monomial basis, the number of prolongations con- 
sidered during completion, the portion of reducible prolongations, and 
the computation time. Thus, 100% reducible prolongations means that 
the input is already an involutive basis. An empty entry in the column 
for Pommaret division means that we did not compute a minimal Pom- 
maret basis because the ideal is not zero dimensional. For the other 
divisions, it means that the timing is larger than 10000 seconds at our 
computer[|. 
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16 


6 
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For some examples, bases for two different divisions may coincide. For 
the fourth cyclic roots (Example |l0|), the bases for Thomas division 
and Division I, as well as those for Janet division and the induced 
division Djjl coincide, respectively. 

The computations with monomial sets should give at least some 
hint to the performance of different divisions in the polynomial and 



differential cases. Prom our experience, Janet division, generally, and 
Induced divisions, sometimes, seem to be the most promising in terms 
of prolongations that have to be considered. Pommaret division - even 
though it is not noetherian - deserves further investigation, because it 
is globally defined and rather "compact", too. 

3 Conclusion 

In addition to the above described implementation of different invo- 
lutive divisions in Mathematica we implemented the completion al- 
gorithm for Janet division (Example |j) in C. In this case an input 
monomial set is represented as an array of lexicographically ordered 
multiindices and its completion is done with respect to the same order. 
This choice of completion ordering was motivated by the monotonicity 
of Janet division with respect to the lexicographical order. The partial 
involutivity of an intermediate monomial set is preserved in the course 
of completion and the time for recomputation of the separations is 
minimized Q. 

The set of nonmultiplicative prolongations to be treated is also 
represented as a lexicographically sorted array of multiindices, that 
provides the simplest way to choose a minimal prolongation. Every 
time an irreducible nonmultiplicative prolongation occurs it is inserted 
in the intermediate monomial set and its nonmultiplicative prolonga- 
tions are inserted in the prolongation set. The determination of their 
position in the sorted arrays is performed using the binary search algo- 
rithm. In so doing, the check of Janet reducibility of the prolongation 
under consideration is done in the course of the position determina- 
tion. This is a rather straightforward procedure that makes use of the 
partition into prefix-groups as defined in Example ^. 

The C implementation was done in GNU C/C++ version 2.81 on a 
100 MHz Pentium computer running Windows 95. The running times 
for examples in the above table are less than 0.01 seconds, except 
Example ||] for n = 4 which took about 5 seconds. 

We plan to extend both Mathematica and C codes to polynomial 
and then to linear differential systems. Whereas the highly flexible 
and easily extensible Mathematica code allows one to experiment with 
different involutive divisions, in the further development of the C code 
we are going to restrict ourselves to Janet, Pommaret and may be 
Induced divisions which are more preferable from the computational 
efficiency point of view. 
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